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Abstract. We consider the Klein-Gordon equation on a star-shaped network composed of n 
half-axes connected at their origins. We add a potential which is constant but different on each 
branch. The corresponding spatial operator is self-adjoint and we state explicit expressions 
for its resolvent and its resolution of the identity in terms of generalized eigenfunctions. This 
leads to a generalized Fourier type inversion formula in terms of an expansion in generalized 
eigenfunctions. 

The characteristics of the problem are marked by the non-manifold character of the star- 
shaped domain. Therefore the approach via the Sturm-Liouville theory for systems is not 
well-suited. 



1. Introduction 

This paper is motivated by the attempt to study the local behavior of waves near a node in 
a network of one-dimensional media having different dispersion properties. This leads to the 
study of a star-shaped network with semi-infinite branches. Recent results in experimental 
physics pH [19] , theoretical physics [H] and functional analysis [8[ [13] describe new phenomena 
created in this situation by the dynamics of the tunnel effect: the delayed reflection and advanced 
transmission near nodes issuing two branches. It is of major importance for the comprehension 
of the vibrations of networks to understand these phenomena near ramification nodes i.e. nodes 
with at least 3 branches. The associated spectral theory induces a considerable complexity (as 
compared with the case of two branches) which is unraveled in the present paper. 

The dynamical problem can be described as follows: 
Let Ni,...,N n be n disjoint copies of (0,+oo) (n > 2). Consider numbers a k ,c k satisfying 
< Cfc, for k = 1, . . . , n and < a\ < a2 < ■ ■ ■ < a n < +oo. Find a vector (ui, . . . ,u n ) of 
functions u k '■ [0, +oo) x N k — > C satisfying the Klein-Gordon equations 

[df - c k dl + a k ]u k (t, x) = , k = 1, . . . , n, 

on N\ , . . . , N n coupled at zero by usual Kirchhoff conditions and complemented with initial 
conditions for the functions u k and their derivatives. 

Reformulating this as an abstract Cauchy problem, one is confronted with the self-adjoint 
operator A = {—c k ■ d% + a k )k=i ... n i n L 2 (N), with a domain that incorporates the Kirchhoff 
transmission conditions at zero. For an exact definition of A, we refer to Section [2j 

Invoking functional calculus for this operator, the solution can be given in terms of 

e^^uo and 
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The refined study of transient phenomena thus requires concrete formulae for the spectral rep- 
resentation of A. The seemingly straightforward idea to view this task as a Sturm-Liouville 
problem for a system (following [27J) is not well-suited, because the resulting expansion for- 
mulae do not take into account the non-manifold character of the star-shaped domain. The 
ansatz used in [27] inhibits the exclusive use of generalized eigenfunctions satisfying the Kirch- 
hoff conditions. This is proved in Theorem 18.11 in the appendix of this paper, which furnishes 
the comparison of the two approaches. 

A first attempt to use well-suited generalized eigenfunctions in the ramified case but without 
tunnel effect [5] leads to a transformation whose inverse formula is different on each branch. 
The desired results for two branches but with tunnel effect are implicitly included in [27] . For n 
branches but with the same Ck and on all branches a variant of the above problem has been 
treated in [7] using Laplace transform in t. 

In the present paper we start by following the lines of [5]. In Section we define n families 
of generalized eigenfunctions of A, i.e. formal solutions for AG [a\, +oo) of the equation 

AF% = AF A fc 

satisfying the Kirchhoff conditions in zero, such that e ±iv ^*i^ (x) represent incoming or outgoing 
plane waves on all branches except for A G [a n , +oo). For A G [a p , a p+ i), 1 < p < n we have 
no propagation but exponential decay in n — p branches: this expresses what we call the multiple 
tunnel effect, which is new with respect to [5]. Using variation of constants, we derive a formula 
for the kernel of the resolvent of A in terms of the Ft . 

Following the classical procedure, in Section H] we derive a limiting absorption principle for 
A, and then we insert A in Stone's formula to obtain a representation of the resolution of the 
identity of A in terms of the generalized eigenfunctions. 

The aim of the paper, attained in Section [7J is the analysis of the Fourier type transformation 



in view of constructing its inverse. We show that it diagonalizes A and determine a metric 
setting in which it is an isometry. This permits to express regularity and compatibility of / in 
terms of decay of Vf. 

Following [3] up to the end would induce a major defect in the last step of this program: the 
Plancherel type formula would read 



where the indices j,j + 1 are considered modulo n, and Kj is a suitable weight. This cyclic 
structure stems from the underlying formula for the resolvent derived in Section [31 which reflects 
the invariance of a star-shaped network with respect to cyclic permutation of the indices of the 
branches and thus the non-manifold character of the domain. This feature inhibits the analysis 
of the decay properties of the (Vf)k'- the finiteness of \\f\\ 2 H does not automatically imply the 
decay of the terms on the right-hand side. In fact, the cutoff by the characteristic function 1-^- 
causes a poor decay in A. 

Consequently, the main objective of the present paper is the elimination of the cyclic structure, 
which is inevitable in the kernel of the resolvent, from the Plancherel type formula. To this end, 
we use a symmetrization procedure leading to a true formula of Plancherel type 
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This is carried out in Section [5] combining the expression for the resolution of the identity E(a,b) 
found in Section [J] with an ansatz for an expansion in generalized eigenfunctions: 



This creates an (3n 2 + 1) x n 2 linear system for the qi m , whose solution leads to the result in 
Theorem 15.31 and to the Plancherel type formula. 

A direct approach to the same symmetrization problem, carried out in Section El yields a 
closed formula for the matrix q based on n x n matrices. This approach is to a great extent 
independent of the special setting and is thus supposed to be generalizable. 

In Section [7J the desired inversion formula as well as the Plancherel type theorem are stated. 
Finally the domains of the powers of A are characterized using the decay properties of Vf. We 
see that V exhibits all features of an ordered spectral representation (see Definition XII. 3. 15, p. 
1216 of [16]) except for the surjectivity, which is not essential for applications. The spectrum 
has n layers and it is p-fold on the frequency band [a p , a p +i). This reflects a kind of continuous 
Zeemann effect due to the perturbation caused by the constant, semi infinite potentials on the 
Nj given by the terms ctjUj. On this frequency band the generalized eigenfunctions have an 
exponential decay on n — p branches, expressing the multiple tunnel effect. 

Our results are designed to serve as tools in some pertinent applications concerning the 
dynamics of the tunnel effect at ramification nodes. In particular, we think of retarded reflection 
(following [HI H]), advanced transmission qt barriers (following |17[ Q31 03]), L°°-time decay 
(following [21 [3]), the study of more general networks of wave guides (for example microwave 
networks [25]) and causality and global existence for nonlinear hyperbolic equations (following 



Finally, let us comment on some related results. The existing general literature on expansions 
in generalized eigenfunctions ( [11} [2^1 [27] for example) does not seem to be helpful for our kind 
of problem: their constructions start from an abstractly given spectral representation. But in 
concrete cases you do not have an explicit formula for it at the beginning. 

In [10] the relation of the eigenvalues of the Laplacian in an L°°-setting on infinite, locally finite 
networks to the adjacency operator of the network is studied. The question of the completeness 
of the corresponding eigenfunctions, viewed as generalized eigenfunctions in an L 2 -setting, could 
be asked. 

In [21j, the authors consider general networks with semi-infinite ends. They give a construction 
to compute some generalized eigenfunctions from the coefficients of the transmission conditions 
(scattering matrix) . The eigenvalues of the associated Laplacian are the poles of the scattering 
matrix and their asymptotic behaviour is studied. But no attempt is made to construct explicit 
inversion formulas for a given family of generalized eigenfunctions. 

Spectral theory for the Laplacian on finite networks has been studied since the 1980ies for 
example by J. P. Roth, J. v. Below, S. Nicaise, F. Ali Mehmeti. A list of references can be found 



In [22j the transport operator is considered on finite networks. The connection between the 
spectrum of the adjacency matrix of the network and the (discrete) spectrum of the transport 
operator is established. A generalization to infinite networks is contained in [15] . 

Gaussian estimates for heat equations on networks have been proved in [23J. 

For surveys on results on networks and multistructures, cf. [H I18j. 

Many results have been obtained in spectral theory for elliptic operators on various types of 
unbounded domains for example [201 HI El E] , cf. especially the references mentioned in [3]. 
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2. Data and functional analytic framework 

Let us introduce some notation which will be used throughout the rest of the paper. 

Domain and functions. Let Ni,... ,N n be n disjoint sets identified with (0, +00) (n G N, 
n > 2) and put N := Ufc=i Nk, identifying the endpoints 0, see [1J for a detailed definition. 
Furthermore, we write [a, b]isr k for the interval [a, b] in the branch Nk- For the notation of 
functions two viewpoints are used: 

• functions / on the object N and fj~ is the restriction of / to Nk- 

• n-tuples of functions on the branches Nk', then sometimes we write / = . . . , f n ). 

Transmission conditions. 

n 

(T ): (uk)k=l,...,n G JJ C °(Nk) satisfies Uj(0) = tifc(0) V i, k G {1, . . . , n}. 
fe=l 

This condition in particular implies that (y,k)k=l,... n may be viewed as a well-defined function 
on N. 

n n 

(Ti): (u k )k=l,...,n G JJ C 1 ^) satisfies J^c fc • ^(0+) = 0. 
fc=i fc=i 

Definition of the operator. Define the real Hilbert space H = Y\j t= iL 2 (Nk) with scalar 
product 

n 

{u,v) H = ^2{uk,v k ) L 2( Nk) 

k=l 

and the operator A : D(A) — ► H by 

n 

D(A) = |(ujfe)fc=i,... in G Y[ H2 ( N k) '■ i u k)k=i,...,n satisfies (T ) and (Ti)\, 

k=l 

A(( u k)k=l,...,n) = (AkU k )k=l,...,n = (~Cfc ■ d^Uk + a k U k )k=l,...,n- 

Note that, if Ck = 1 and = for every k G {1, . . . , n}, A is the Laplacian in the sense of the 
existing literature, cf. flO ] I2T1 123] . 

Notation for the resolvent. The resolvent of an operator T is denoted by R, i.e. R(z,T) = 
(zI-T)- 1 for z G p(T). 

Proposition 2.1. The operator A : D(A) — ► defined above is self-adjoint and satisfies o~(A) C 
[ai, +cxd). 

Proof. Consider the Hilbert space 

V = < (uk)k=l,...,n G H l (N k ) : (wjfc)fe=i ) .„,„ satisfies (T ) ^ 
^ fe=l > 

with the canonical scalar product (•, -)y. Then the bilinear form associated with A + (e — a\)I 
is a £ : V x V — > C with 

n 

a e (u,v) = y~] [c fc (g x itfc, d x Vk)L2 (N k ) + (a fc + e - ai)(u fc , ■u fc ) L 2( A r fc )] . 
fe=l 

Then clearly there is a C > with a e (u,u) > C(u,u)y for all u G ^ and all e > 0. By partial 
integration one shows that the Friedrichs extension of (o e , V, H) is (A + (e — a{)I, D(A)). Thus 
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the operator A + (e — a\)I is self-adjoint and positive. Hence a(A + (e — a\)I) C [0, +00) for all 
e > 0, what implies the assertion on the spectrum. □ 

3. Expansion in generalized eigenfunctions 

The aim of this section is to find an explicit expression for the kernel of the resolvent of the 
operator A on the star-shaped network defined in the previous section. 

Definition 3.1. An element / G 11^=1 C°°(Nk) is called generalized eigenfunction of A, if it 
satisfies (Tq), (Ti) and formally the differential equation Af = Xf for some A G C. 

Lemma 3.2 (Green's formula on the star-shaped network). Denote by i n the subset of the 
network N defined by 

n 

V h ,...,i n = {0}u{J(0,l k ) Nk . 

k=l 

Then u, v G D(A) implies 



n n 

u"(x)v(x) dx = I u(x)v"(x) dx — u(lk)v'(lk) + u'(lk)v 

(_ JVl, l„ 1 1 ; 1 



Proof. Two successive integrations by parts are used and since both u and v belong to D(A), 
they both satisfy the transmission conditions (To) and (Ti). So 

n n 

5>*(0K(0)=ui(0)J>U0) = 0. 

k=l k=l 

Idem for ELi<(°K( )- D 

This Green formula yields now as usual an expression for the resolvent of A in terms of the 
generalized eigenfunctions. 

Proposition 3.3. Let A G C with Im(A) ^ be fixed and let e\, e\ be generalized eigenfunctions 
of A, such that the Wronskian Wi 2 ( x ) satisfies for every x in N 

wl 2 (x) = det W{e\{x), e x 2 (x)) = e\{x) ■ (e x 2 )'{x) - {e\)'{x) ■ e x 2 (x) + 0. 

If for some fe € {1, ... ,n} we have e 2 \N k G H 2 (Nk) and e^|Ar m G H 2 (N m ) for all m ^ k, then 
for every f G H and x G Nk 



[R(X,A)f](x ' 



c k w^ 2 (x) 



x (x)e x (x')f(x') dx' + f e x {x)e x {x')f{x') dx' 

(x,+oo) Nk J N\(x,+oo) Nk 



(1) 



Note that by the integral over N, we mean the sum of the integrals over Nk, k = 1, . . . , n. 



Proof. Let A G p(A). We shall show that the integral operator defined by the right-hand side of 
([T]) is a left inverse of XI — A. Let u G D(A) and x G Nk- Then 

I\ := J e x (x)e x (x')(XI - A)u{x') dx' + f e x (x)e x (x')(XI - A)u{x') dx' 

J(x,+oo) Nk JN\(x,+oo) Nk 

= e x {x) lim [ " e 2 (x')(XI - A)u(x') dx' 

+ el{x) lim / e\{x'){XI - A)u{x') dx' , 

l m ^oo,m^k Jy 
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We have u G D(A) C U]=i H 2 ( N j) and 



due to the dominated convergence Theorem, the integrands being in L (R) by the hypotheses. 

n-=i^ 2 (^)and 

e$\ Nk G H 2 (N k ), e\\ Nm G H 2 (N m ), m^k 



by hypothesis and thus 

d x u\ Nk (x) ■ e 2 \ Nk (x) — > 0, u\ Nk (x) ■ d x e 2 \ Nk (x) — ► 0, 

x — >+oo " x—>-\-oc 

d x u\ Nm {x) -e\\ Nm {x) —> 0, u\ Nm {x) ■ d x e\\ Nm (x) — > 0,m^k, 

x— >+oo x^+oo 

all products being in some H 2 (Nj). Recall that 

' fg = I fg" - f(b) 9 '(b) + f'(b)g(b) + f(a)g'(a) - f'(a)g(a) 



for f,g e H 2 ((a, b)). Now Lemma [3721 and (A/ - A)e$ = for r = 1, 2 imply 
/ A = ei(x) lim \ f * (XI - A)e%(x')u(x') dx' 



+ c k [~u(l k )d x e 2 (l k ) + d x u(l k )e 2 (l k ) + u(x)d x e 2 (x) - d x u(x)e 2 (x) 
+ e%(x)\ lim / (XI - A)e\{x')u{x) dx 

ll m ^oo,m^k Jy 

+ ^2cj lim (-u(lj)d x e^(lj) + d x u(lj)ei(lj)\ + c k ( y -u(x)d x e\ (x) + d x u(x)e\ (x) ^ 

= c k (e\(x)d x e 2 (x) - d x e\{x)e 2 {x) S ju(x) 
= c k Wi 2 (x)u(x). 

Now the invertibility of XI — A implies the result. □ 
Definition 3.4 (Generalized eigenfunctions of ^4). For k G {1, . . . ,n} and A G C let 



c / /A-a fc Y,i^k c &( X ) 

^ (A):= V^ Sfc:= " c fc &(A) ■ 

Here, and in all what follows, the complex square root is chosen in such a way that V r • e 1 ^ = 
^/f e ^/2 w it]a r > and G [— 7T, 7r). 

For A G C and j, k G {1, . . . , n}, F^' j : N -> C is defined for x G iV^ by T^' j (x) := F^ k j (x) 
with 

f /-.v/i.r) = coste(A)z) ± ^(A)sin(^(A)x), 

I F x,k( x ) = exp(±i^(A)x), for fc ^ j. 

Remark 3.5. • i*\ ^' satisfies the transmission conditions (To) and (Ti). 

• Formally it holds AF^' j = XF^' j . 

• Clearly T^'- 5 does not belong to H, thus it is not a classical eigenfunction. 

• For Im(A) ^ 0, the function F^'J, where the +-sign (respectively — sign) is chosen if 
Im(A) > (respectively Im(A) < 0), belongs to H 2 (N k ) for k ^ j. This feature is used 
in the formula for the resolvent of A. 
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Figure 1. NxN in the case n = 3 



Definition 3.6 (Kernel of the resolvent). Let w : C — > C be defined by w(X) = ±i'X)^=i c j£j(^)- 
For any A G C such that w(A) ^ 0, j 6 {1, . . . , n} and x £ Nj, we define 

— ^-yF^(x)F^ +1 (x'), for x' G iV,-, x' > x, 

7v\F\j +1 ( x )F\ ' J \ x ')i f° r x' E Nk,k ^ j or x' E Nj, x < x. 
(A) 

In the whole formula + (respectively — ) is chosen, if Im(A) > (respectively Im(A) < 0). Finally, 
the index j is to be understood modulo n, that is to say, if j = n, then j + 1 = 1. 

Figure [T] shows the domain of the kernel K(-,-,X) in the case n = 3 with its three main 
diagonals, where the kernel is not smooth . Note that in particular, if Cj = c and aj = for all 
j G {1, . . . , n}, then w(X) = ±mc\/A for all j G {1, . . . , n}, which only vanishes for A = 0. On 
the other hand, if there exist i and j in {1, . . . , n}, such that aj 7^ aj, then it is clear that w(X) 
never vanishes on R, but we need to know, if it vanishes on C. 

We will show in Theorem 13.81 that K is indeed the kernel of the resolvent of A. In order to 
do so, we collect some useful observations in the following lemma. 

Lemma 3.7. i) For a\ < X and e > holds \w(X — ie)\ 2 > X^j=i c j\^~ a j\- 

ii) For X G p(A) such that Re(A) > a\, the Wronskian w only vanishes at A = a, if a& = a 
for all k G {1, . . . , n}. 

Proof. We first prove H]). 



K(x, x , A) = < 
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Note that for z\ , zi , . . . , z n G C holds 



3=1 3=1 



k,l=l 
k^tl 



With z k := c k ^ k (X — ie) and the abbreviation rj k := y/X — ie — a k it follows 



(A - ie)\ 2 =^2 Cj \rjj\ 2 + 2 ^ c k ciRe(r] k r]i) 

3=1 



k,l=l 
k^l 



Thus it suffices to show Re(??fc7/;) > for fc, Z = 1, . . . , n with k ^ I. With our convention 
^/i = y/\z\e l ~^~ , arg(z) G [— 7r, 7r), this means 



arg(%r/i) G 



7T 7T 

2' 2 



Without loss of generality, let k < I. Then we have a k < a/ and there are three possible positions 
of A: 

• ai < A: 

Then for r = fc and for r = Z we have arg(A — ie — a r ) G [—§,0] and therefore 



arg(r?r) = ^ ar s( A ~ ie ~ a r) G 
Using A — afc > A — a/, this implies 

arg(r) k rji) = arg(r] k ) - arg(r#) G 

• a k < \ < (if. 

Then arg(A — ie — a k ) G [— f ,0] and therefore 

arg(%) = - arg(A - ie - a k ) G 



7T 

-4'° 



7T 



7T 



-,o 



Furthermore, arg(A — ie — ai) G [—it, — |], and therefore 



arg(?7/) = - arg(A - ie - a t ) G 



7T 7T 

2'~4 



Putting everything together we get 



arg(%7ft) = arg(r/ fe ) - arg(^) G 



0, 



7T 



• A < a k : 

In this case we get again for r = k and r = I arg(A — ie — a r ) G [—it, — |] and therefore 



This yields 



arg(r? r ) = - arg(A - ie - a r ) G 



^g(iJkVl) = arg(%) - arg(rft) G 



7T 7T 

2'~4 

7T 7T 

~4' 4 



Thus, in all three cases we have arg^r^) G [—§,§] and, hence, Re(r) k r]i) > for all fc,Z = 
1, . . . , n with k 7^ Z. 
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In order to prove [ITJ) , we note that the choice of the branch cut of the complex square root 
has been made in such a way that = y/Jl for all /J, G C. Thus w(jl) = w(p) for all fi G C. 
This implies that the first part of the lemma can be generalized to: 

n 

\w(p)\ 2 > ^Cj|Re(» - aj\ 

for every \i such that Re(/x) > a\. Then w(fi) = 44> |Re(/i) — dj\ = for all j G {1, . . . , n}. 
This implies that to never vanishes, if there exist k, I in {1, . . . , n}, with k ^ I such that a k ^ d[. 

Finally, if au = a for all k G {l,...,n}, then ^(/z)! 2 = (J2k=iV^k) \p — ct\, which only 
vanishes for p = a. □ 

Theorem 3.8. Let f G H. Then, for x G N and A G p{A) such that Re(A) > a 1 

[R(X,A)f](x)= [ K(x,x',X)f(x') dx'. 

JN 

Proof. In ((T]), the generalized eigenfunction e\ can be chosen to be i*\ Then can be -P 1 ^'' 
with any I ^ j so we have chosen j + 1 to fix the formula. The choice has been done so that 
the integrands lie in L 1 (0, +oo) (cf. the last item in Remark 13. 5|) . Simple calculations yield the 
expression for the Wronskian u;(A). 

For A G p(A) such that Re(A) > ai, the Wronskian only vanishes at A = a if = a for all 
k G {1, . . . , n} due to Lemma [3771 But in this case, the Wronskian is w(X) = (T2k=l V^k) VX — a 
and w(X)~ 1 has an L 1 -singularity at A = a. □ 

4. Application of Stone's formula and limiting absorption principle 
Let us first recall Stone's formula (see Theorem XII. 2. 11 in |16]). 

Theorem 4.1. Let E be the resolution of the identity of a linear unbounded self-adjoint operator 
T : D(T) — > H in a Hilbert space H (i.e. E(a,b) = l( a ^(T) for (a,b) G M?, a < b). Then, in 
the strong operator topology 

h(T)E(a, b) = lim lim — /" h(X)[R(X - ei, T) - R(X + ei, T)} dX 

5^0+ e^0+ 2TTI J a+5 

for all (a, b) G M 2 with a < b and for any continuous scalar function h defined on the real line. 

To apply this formula we need to study the behaviour of the resolvent R(X, A) for A approach- 
ing the spectrum of A. Lemma 13.71 will be useful as well as the following results. 

Lemma 4.2. Let 5 > be fixed. For all a\ < A, < e < <5 and j = 1, . . . , n holds 

1 n 

\s 3 (X -ie)\< M(A, 5) := . max { _ } £((A - a k f + 5 2 f'\ 

Proof. This follows directly from the definition of Sj. □ 
Note that M(-, 5) G L 1 1 oc ([ai, +oo)). Furthermore, if a\ = . . . = a n , then 

I n 
Sj{p) = ^=X)Vcfc 
V c i fe =i 
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for all fA G C, which means that Sj is constant and we may take 



1 n 

M{\,8):= max {— VV^j. 



■3 k=l 

Theorem 4.3 (Limiting absorption principle for A). Let 5 > be fixed and let M(X,5) be 
defined as in Lemma \4-2\ Then for all a\ < X, < e < 5 and (x,x') G N 2 we have 

i) lim a ^o K(x, x' , A — ia) = K(x, x', X), 

ii) \K(x,x',X-ie)\ < N(X, 5)e^ x+x '\ where N(X,5) : = -^i±M^)_ ana > 

(2^j=i c jl A a j\) 

_ i x 
7 := max J= i i ... jn {c i 2 }max{((a n - a x ) 2 + 5 2 )4, l,<5}. 

Proof. i) The complex square root is, by definition, continuous on {z G C : Im(z) < 0} (cf. 

Definition I3.4D . hence the continuity of K(x, x' , •) from below on the real axis. Note that 
x and x' are fixed parameters in this context, 
ii) For Im(/i) < 0, /j, = X — is and x G iVj we have in concrete terms 

[cos(£j(/i)x) — isj(n) sin(£j(//)x)] exp(— i^-(^)x'), x' G iV,, x' > x, 

< exp(— i^j(n)x) [cos(£j(/x)x') — isj{n) sin(^(^)x')] , x' G Nj, x < x, 

k exp(-i^(//)x) exp(-^ fc (/i)x'), x' G fc / j. 

Now, let us first look at the case A > a n . Then 

\exp(-i£j((j,)x)\ < exp(|Im(^-(/x)x)|) = exp(c~ 1/2 |Im( y'A - ie - aj)\). 

Using the fact, that for z G C with | arg(z)| < ir/2 we have |Im(y^)| < max{l, |Im(z)|}, 
we obtain 

| exp(— i£j(p)x)\ < exp(c J - max{l, <5}x) . 
In the case a\ < X < a n we find 



K(x, x 1 , fj,) 



w(fi) 



M»)\ = J lj ^ = Cj - 1/2 ((A - aj f + e 2 )^ < c^\(a n - a,) 2 + S 2 f\ 



Using these estimates and \e z \ = e Re ^ < e' 2 ', we find for all A > a\ 

1 J i 1 + M^)l) exp(|^(^)|x)exp(|^(/i)|x'), x' G Nj, 
\w(n)\ \ exp(|^(/i)|x) exp(|4(^)|x'), x' el%,k^j 

1 



\K(x,x',n)\ < 



< 



[1 + exp( 7 (x + x')). 



The conclusion now follows using Lemma 13.71 and Lemma 14.21 □ 

Note that these estimates in particular imply that N(-,5) G Ll oc ([ai, +oo)). In fact, if a\ = 
. . . = a n , then M(X,8) can be chosen to be constant, see Lemma 14.21 an d the denominator of 
N causes only an Lj oc -type singularity. On the other hand, if there are two different aj, the 
denominator of N is never zero and M(-, 5) G Lf ([ai, +oo)) again by Lemma 14.21 

Lemma 4.4. For (x,x') G N 2 and X G C, it holds K(x,x',X) = K(x,x',X). 

Proof. The choice of the branch cut of the complex square root has been made such that VA = 
VX for all A G C. This implies e^ x = e iy/lx = e '^ x for all A G C and i£l, Thus it holds 

F+' j (x) = F-T J (x) and F~' j (x) = F±' j (x) 
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for all A E C, x E iV and j E {1, . . . ,n}. In the same way we have w(X) = —w(\). Observe, 
that switching from A to A the sign of the imaginary part is changing, so in the definition of 
K(x, x' , A) we have to take the other sign, whenever there is a ±-sign in the formula. This gives 
the assertion. □ 

Now, we can deduce a first formula for the resolution of the identity of A. 

Proposition 4.5. Take f E H = YYj=i L 2 (Nj), vanishing almost everywhere outside a compact 
set B C N and let — oo < a < b < +oo. Then for any continuous scalar function h defined on 
the real line and for all x £ N 

(h(A)E(a,b)f)(x) 



= [ K\) £ / /OO 00 • lm(-^F-^x)F-^ 1 (x')) 

J (a,&)n[ai,+oo) JT[ 3 <J L 3 \W(A) ) 

+ l^+oo)^ 00 ■ ^{^F^ +1 (x)F^(x'))] dx'} dX, 

where E is the resolution of the identity of A (cf. Theorem and the index j is to be 
understood modulo n, that is to say, if j = n, then j + 1 = 1. 

Proof. The proof is analogous to that of Lemma 3.13 in [2]. Let g E H be vanishing outside B. 
Then 

rb-8 



lim 

8,e^o+ 2-Ki 



(h(A)E(a, b)f, g) H =( lim lim / h(X) [R(X - ei, A) - R(X + ei, A)] dX /, g] (2) 
\<5->o+ e^o+ 2m J a+S J h 

lim lim — f ( h(X)\R(X-ei,A)- R(X + si,A)] d\f,g] (3) 
5 ^o+ £ ^o+ 2m \J a+s J h 

1 f b ~ 5 

lim lim / h(X)(\R(X-si,A)-R(X + ei,A)]f,g)„dX (4) 

5^0+ e^o+ 2m J a+5 

hWU f(x')[K(-,x',X-ie) - K(-,x',X + ie)] dx',g(-)\ dX 

H (5) 

KX)(J f(x') [K(-,x, X - is) - K(-,x', X - ie)] dx', g(-)\ dX 

H (6) 

h(X)( [ f(x')2ilm(K(.,x',X-ie))dx',g(-)) dX (7) 

8,e-*0+ lm J( a +S,b-S)n[a u +oo) \J N Jh 

lim - [ h(X)( [ f{x')[ lim lm(K (■ , x' , X - ie))] dx',g{-) j dX (8) 

5^0+ VT J (a+5,b-8)n[ai,+oo) \J N e^0+ J H 



(a+8,b-8)n[ai ,+oo) 



lim 

8,e^o+ 2m 



lim 



(a+8,b-8)n[ai ,+oo) 



(a,6)n[oi,+oo) JN 



h{X) f(x')lm(K(-,x',X-i0)) dx' dX,g(-)\ (9) 



H 



[If h(X){[ /(^Imf-^-^ilv-Wf 1 ^-) N Xx')F^{x)F-^\x') (10) 

JN K ^(a,6)n[ai,+oo) K JN lw \^) ~[ 3 V 3 

+ lN\(x,+oo) Nj ( x ') F \' 3+1 ( x ) F \ ,J ( x/ ))] dx'} dXg(x) dx. 



Here, the justifications for the equalities are the following: 
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©I Stone's formula (Theorem I4.ip . 

(J3J) : After applying the operator valued integral to /, the two limits are in H. So they 
commute with the scalar product in H. 

("/)5)h * s a continuous linear form on C(H), and can therefore be commuted with 
the vector-valued integration. Note that A i— ► R(X, A) is continuous on the half-plane 
{A G C : Re(A) < a\}, since the resolvent is holomorphic outside the spectrum, cf. 
Proposition 12.11 
([5]): Theorem | 
Lemma [ 

z — ~z = 2i ■ Im z \/z € C. 

Dominated convergence. Since supp /, supp g and [a, b] are compact, we use the limiting 
absorption principle (Theorem 14.30 . 
(|9|): Fubini's Theorem. 
(fT0|) : Definition E! 

The assertion follows, since g was arbitrary with compact support. □ 

The unpleasant point about the formula in the above proposition is the apparent cut along 
the diagonal {x = x'} expressed by the characteristic functions in the variable x' . In fact, there 
is no discontinuity and the two integrals recombine with respect to x' . This is a consequence of 
the next lemma that gives an explicit representation of the integrand above. 

In the following, we use the convention 

a n +l ■= +oo, 

in order to unify notation and we set 

S : 'V 

Lemma 4.6. Let j,k £ {1, . . . , n} and let A be fixed in (a p ,a p+ i), with p £ {1, . . . ,n}. Then 
Im[^(i^' J+1 ) .(x)(i^' J ) fc (ic')] is given by the following expressions, respectively: 

• Ifk>j>porj>k>p (Case (a)) 

ImflV^-^', 

• Ifj<k<pork<j<p (Case (b)) 

Im^— ^ cos(£j-x) cos^x') — Im(— ^ sin(£jx) sin(^x') 

— Re^— ^ cos(£j-x) sm(^kx') — Re^— ^ sin(^x) cos(^fcx'), 

• Ifj = k<p (Case (b), j = k) 

Im^— ^ cos(^jx) cos(^x / ) — Im^— ^ sin(^x) sin(^x') 

— Ref — \ cos(£ 7 x) sin(^fcx') — Ref — \ sin(£ 7 x) cos(^x'), 

• If j <p<k (Case (c)) 

Im(— ^je~^ x ' cos(^x) + Im^— je"^ 1 ' sin(^x), 

• Ifk<p<j (Case (d)) 



lm(—\e ^ x cos(^fcx') — Imf — )e ^ z sin(^x'). 
V w / V iw J 
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Proof. Since A belongs to (a p ,a p+ i), £j(A) is a real number, if and only if j < p. Otherwise, £j 
is a purely imaginary number and we have £j =: — 

Now, the proof is pure calculation, using the following expressions for the generalized eigen- 
functions in the case j 7^ k 



lm{F-' j ),(x) 



if £fc G iR, i.e. k > p, 

-suifax), if £ fc G R, i.e. k < p, 
0, if £,k G iR, i.e. fc > p, 



and for j = k 

R fp-'h ( \ -l cos (£kx) +Im.(s k )sm(£ k x), if £ k G R, i.e. fc < p, 

ne l^A J fc W \ Re (i(i + Sfc )) e -«k* + Re (|(1 - s k )) e*** if ^ G iR, i.e. A; > p, 

T (F~' k \ ( \ = f _Re ( s fc) sin (^ x )' if£fcGR, i.e. k < p, 

Lm \*\ ) k \ x ) \ i m (i(i + Sfe )) e -ftx + Im (i(i _ Sfe )) if £ fc 6 iR) j e fc > pj 

respectively. □ 

Theorem 4.7. Ta/ce f € H = YYj=i L 2 {Nj), vanishing almost everywhere outside a compact set 
B C N and let —00 < a < b < +00. Then for any continuous scalar function h defined on the 
real line and for all x G N 

(h(A)E(a,b)f)(x) = [ h(X)J2^){f f(^\^F^ + \^Fx d ^)]dx'}d\, 

(11) 

where again the index j is to be understood modulo n, i.e, if j = n, then j + 1 = 1. 

Proof. All the work has already been done. It remains to inspect the formulae for j = k in the 
cases (a) and (b) of Lemma 14.61 to observe that the expressions are symmetric in x and x' . So 
for x,x' G Nj with x < x' we hnd F7^ (x)F7^ + [x') = F x '• 3 ] (x')Fl (a;), which implies the 
assertion. □ 

5. Symmetrization 

As was already explained in the introduction, the aim of this section will be to find complex 
numbers qi t m, l,m £ {1, . . . , n}, such that the resolution of identity of A can be written as 

(E(a,b)f) (x) = f qim(X)F^ l (x) [ F^(x')f(x') dx' dX, (12) 

Ja Lm=l jN 



in order to eliminate the cyclic structure of the formula in Proposition 14.51 

In this section we shall often suppress the dependence on A of several quantities for the ease 
of notation, so Sj = Sj(X), q l>m = qi, m (X), = ^(A), Q = ^(A) and w = w(X). 

Lemma 5.1. Given x G Nj, equation ()12p is satisfied for all a\ < a < b < +00 and all 
f G L 2 (N) with compact support, if and only if for all j = 1, . . . , n 



Im 



l-F-^ 1 (x)F-^x')] = <T q lm (X)F-\x)F-' m {x') (13) 

l,m=l 

for almost all x 1 G N and X > ai . Here again the index j has to be understood modulo n, that 
is to say, if j = n, j + 1 = 1. 
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Proof. If functions qi jTn , I, m = 1, . . . , n, satisfy (|T3|) and if a\ < a < b < +00, we get by ([IT 

b n n 

' dA 



/It, n I h 

EM»0{ f{x') qi m (X)F-\x)F~^(x')dx'} 

a i =1 N l ' m=1 

= / E Qlm(X)F-'\x) [ F^[x')f{x') dx' dX, 

Ja l,m=l JN 



which is (fT2j) . 

For the converse implication, let (I12p be satisfied for some x E Nj and all a\ < a < b < +00, 
as well as all / E L 2 (N) with compact support. This means by (jlip 



f(x )Im 



A? 



^ +1 (x)F A -"*(xO 



dx' dA 



:,m=i 



E ^(A)F--'(x) / F A - m (x')/(x') dx' dX. 



Firstly, we want to see that the integrands of the A-integrals on both sides in fact have to be 
equal almost everywhere. In order to do so, we observe that they both are in L 1 ((ai, +00)) and 
use the following general observation: If I is an interval and g E L l (I) satisfies fjg = for 
all intervals J C J, then g = almost everywhere in I. Indeed, in this case we have for any 
Lebesgue point xq E / of g and every e > (cf. [26l Theorem 8.8]) 



I^Oo)! 



1 

Ye 



X +£ -y 

(g(xo) - g(x)) dx + 

X -£ 



2s 



xo+e 



g[x) dx 



X -£ 



1 

< 

~ 2s 



X Q +£ 



\g(x) - g(x )\ dx — >0 (e 0), 



X -£ 



which implies g(xo) = for almost all xq E I. 
This implies 



1 



N 



f(x')lm -^-F^ +1 (x)F^(x') 



dx' 



N 



E q lm (X)F-' l (x)F-' m (x')f(x')dx' 



Lm=l 



for almost all A > a\ and all / E L 2 (N) with compact support. By the fundamental theorem of 
calculus this implies the assertion. □ 

In a next step, we explicitely write down equation fjl3[) as a linear system for the values q\ m . 
Lemma 5.2. The equation 
1 



Im 



HF^)ix){F-^ 



w 



E * m (A)(F A -%.(x)(F A -' m ) fc (x') 



i,m=l 



holds for any (x,x') E Nj x and A E (a p ,a p +i), «/ and onh/ if 
• Case (a): if k > j > p or if j > k > p 

qjk = 

= 

2~2m^k Ijm = 
, J2l^j,m^k llm = I m [^] ! 
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Case (b): if j < k < p or if k < j < p 



J2l^j,m^k Qlm + s fe Yl/l^j Ilk 
J2l^j,m^k llm + s k Qlk 
< J2l^j,m^k llm + Yltyj Ilk 

Case (b), j = k: if j = k < p 

J2lytj,m^j Um +Sj llj 
J2l^j,m^j Urn +Sj Uj 
. J2ljLj,m^j llm + llj 

Case (c); if j < p < k 



+ S i Em^fc Ijrr 
+ Em^lc Qjm 
+ S j J2m^k Qjn 



+ S i Em^' Ijr 
+ Em^j Ijm 
+ S j Em^j Ijr 



^2l,m llm 
+Sj ■ Sfc • qjk 

• Qjk 
+Sj ■ q jk 



z2l,m Qlr. 

+ Sj ■ Sj 

+ s i ' Ijj 
+ s i ' Ijj 



Ijj 



I- [£] 
-I- [h] 
~i ■ Im [&] 

-i • Im \±] 
i ■ Im [4-1 , 



qjk 
J2i^j iik 

, J2l^j t mjtk Urn + Sj J2 m ^k 9jm 

Case (d) : if k < p < j 



. ^2l^j,m^kUr. 



qjk 

J2l^j Uk 
Ylm^k 1j m 
+ SkJ2l^j Uk 



-- 
: 

= i • Im r J- 

Ljui 




-* • Im [in 



Proof. The sum Y17m=l QimW{F\ ( X ){F\ m )fc( x ') * s explicitely written in the different cases 
(a), (b), (c) and (d) as it was done in Lemma [4,61 Then the linear independence of the following 
families of functions is used to get the above systems for the qims: 

• Case (a): e Ax+Bx ' , e Ax ~ Bx ' , e ~ Ax+Bx \ e - Ax ~ Bx \ 

• Case (b): cos (Ac) cos(Bx'), cos(Ae) sm(Bx'), sin(Ac) cos(Bx'), sin(Ac) sin^x'), 

• Cases (c) and (d): e Ax cos(Bx'), e Ax sin(Bx'), e~ Ax cos(Bx'), e~ Ax sm(Bx'), 
where A and B are any fixed non- vanishing real numbers. 

On the other hand: if the qi m satisfy the system indicated in the lemma, on both sides of the 
equation of the lemma we have the same linear combination of the functions given above. Thus 
equality holds. □ 

Having the linear systems in the above lemma at hand, it remains to solve them. Indeed this 
is possible and we briefly indicate the necessary steps. 

Firstly, if A < a±, we have £j £ iSL for all j £ {1, . . . , n}. Thus, for all j, k £ {1, . . . , n} and all 

(a?, x') £ Nj x N k 



Im 



1 

-i 

iw 



Fx' j+1 )j^)(^' j ) k (x') 



Im 



-i'x-t'.x' 



o. 



Now, let A be fixed in (a p , a p+ i), with p £ {1, . . . , n}, remembering the convention a n+ \ = +oo. 
Due to the first equation of the corresponding system of cases (a), (c) and (d), the matrix 
q '■= (qim)f m= i has to be of the form 



Qp 
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with a p x p matrix Q p . 

Hence, the equations Y^tyj Qlk = and Ylm^k Qjm = in case (a) are obviously fulfilled, as 
well as YltyjQlk = in case (c) and Ylm^k Qjm = in case (d). This means that only three 
equations remain from the cases (a), (c) and (d): 

£z, m 92m = M£] 

E m ^fc+^Eto = -*-M?5f] 

E^i Qlm + sj £ q jm = Mm [±] , 

where in all the sums, I and m belong to {1, . . . ,p}. Now it is important to note that these 
three equations are already contained in the following system corresponding to the case (b), 
i.e. the only conditions for the qi m 7 s are the four following equations for a fixed (j, k) such that 
j < k < p or k < j < p: 



^l^j^mjkk Qlm 
J2l^j,m^k Qlm 
< ^l^j,m^k Qlm 

and, for j = k < p: 



+~$k Yll^j Qlk 
+~$k Yjtyj Ilk 

+ E^y Qik 



+ s j Em^ Qjr 
+ Em^S; Qjm 
+ s j J2m^k Qjr 



E/ m Qlm 



+Sj ■ Sk ■ qjk 
+~sh~ • qjk 
+sj ■ q jk 



T,l^j,m^j Qlr 
^m,m+j Qlr 
< 12l^j,m^jQlr 



+ s j Em^j Qjr 
+ Em^j Qjm 
+ S i Em/j Qjr 



Ei,m Qlm 
+Sj ■ Sj ■ qjj 

~>~ S j ' Qjj 

+ s i ' Qjj 



+SjI2i^jQij 
+~jHi^jQij 
+ Ej# Qij 

where, once more, / and m belong to {1, . . . ,p} in all the sums 

Now, we denote by Q jk := J2i^jQik and Q' jk := J2 m ^kQjm 
Re \j- 1 , the above system (fT4"l) is 

J2l^j t m^k Qlm 

(Sk ~ l)Qjk+ (Sj - l)Q' jk + (SjSk- - l)q jk 

(sk-l)Q jk + (Sfc-l)gjfc 

( s j - l )Q'jk+ ( s j - l )Qjk 



M£] 

-* • Im lh 

i ■ Im \X\ 

Ma 

-MS] . 

— i • Im T J- 
t • Im [ J-l , 



(14) 



(15) 



Using; the fact that Im [- 



Im [h] ~ Qjk ~ Q'jk - Qjk 



Tin [-= 

L ? 



iw ' 



iw 



Since s,- 



iw 

'/j' 



the last three equations may be rewritten as 



{-ifjw)Q jk + {ifkw)Q' jk + (-ifjW + if k w -\w\ 2 )q jk = (2/j/ fc )Re [^] 



Qjk + 



Q' jk + 



Qjk 
Qjk 



fk 



and the Gauss method yields that c/jfc must vanish for j ^ fe. 

In the case j = k, we rewrite system (I15p as above and the three equations to be solved turn 
out to be 



(-ifjw)Qjj+ (ifjw)Q' jj + (-ifjW + ifjW -\w\ 2 )qjj = {2fjfj) (Re 



Qjj+ 



Q'jj+ 



Qjj 
Qjj 



fi 
h 



The Gauss method once more gives the only possible solution qjj = for any j G {1, . . . ,p}. 

With this candidate for a solution at hand, it is pure calculation to show that it is indeed a 
solution. Thus we have shown the following result. 
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Theorem 5.3. Let A be defined as in Section® and the generalized eigenf unctions F^'° be given 
by Definition \&.4\ Take f G H = n?=i L 2 (Nj), vanishing almost everywhere outside a compact 
set B C N and let —oo < a < b < +00. Then for all x £ N 



j-b n 

(E(a,b)f)(x 



y>(A)F-'<(x) / F-\x')f(x')dx' dX, 
ti Jn 



where 



JO, if\<ai, 
«W := \ fiW if ai< \ 

Furthermore, for almost all X £l the matrix qi tm (X) = <^ m l( aj)+00 )(A)/z(A)/|w;(A)| 2 is the unique 
matrix satisfying 1)120 . 



6. A DIRECT APPROACH 

We have seen in the preceding section that a matrix g(A) satisfying (I12p exists and is unique 
up to a null set. It is the aim of this section to deduce an alternative representation for q(X), 
involving only n x n-matrices and not an (3n 2 + 1) x n 2 system as above. Since this approach is 
essentially independent of the special setting, it should be more convenient for generalisations. 

In the following, let us consider the complex- valued generalized eigenfunctions F7' k as func- 
tions on iV and not as elements of IIj=i C°°(Nj). We introduce the notation 

/ F7' 1 (x) 



Denoting by = (Sik)l=i,...,n = (0, . . . , 0, 1, 0, . . . , 0) T the A:-th unit vector in C n , we set e?i(A) := 
F\(0) and for j = 2, . . . , n we fix Xj G Nj and set dj(X) := F\(xj). Due to the form of our 
generalized eigenfunctions, we then have 

n 

d 1 (A) = ^e fc = (l,...,l) T , 
k=i 

dj(X) = fyej + aij ^ e k = (ay, . . . , aj,/3j,aj, a>j) T for j = 2, . . . , n 



for suitable otj,f3j G C. 

Using these vectors, we now define 



D(X) :=J>(A)e 



T 
3 



Since ay 7^ /3j, for every j G {2, . . . , n}, by construction, the matrix D is invertible for any choice 
of (x2, • • • , x n ) provided that Xj ^ for all j G {2, . . . , n}. Indeed det D = X\^ =2 {fij — ay)- 
Denoting by C(A) the diagonal matrix with cn = i and c,y = ictj for any j G {2, . . . , n}, we can 
formulate our theorem. 
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Theorem 6.1. The matrix q(X) := (qim(^))l,m=i,...,n satisfying (fT2|) is given by 

1 n 

q(X) = (D(Xy 1 ) T lm(^—— J2eAW T ej+iejm)(W) f 

j=i 

= (D(A)- 1 flm(^C'(A) J D(A))(^(A)" 1 ) T 

for almost all A > a\. As previously, j has to be understood modulo n, that is to say, j + 1 = 1 
if j = n. 

Proof. By Lemma 15.11 the function q satisfies (|13p . Using the matrices and vectors introduced 
above, for fixed j G {1, . . . , n} this can be rewritten as 

1 

-w(A) J 

for x G Nj and x' G N. 

Setting x = Xj and x' = x k we obtain by the definition of dj 



Iml—-F x (xYe j+1 ejF x (x')) = F x (xY q(X) F x (x') 



\w.{^-d J {X) T e j+l e]d k {X)) = d j (X) T q(X)d k (X) 



M>) 

for all j, k £ {1, . . . , n} and thus 

n n 1 n n 

k=lj=l ^ ' ' k=lj=l 



Using Ylk=i^k^k = D, we obtain 



i n 

Y,^ d M) T e J+ ie^D) = D T qD 



Im, 

M*) ■ 

Since D is invertible and the relation YHj=i e jdj{X) T ej+ieJ = —iC(X) holds true, we now get 
the desired formula for q. □ 

7. A Plancherel-type Theorem 

In this section we prove a Plancherel type theorem for a Fourier type transformation V and 
its left inverse Z associated with the generalized eigenfunctions (i ? A ^' fc )fc=i,...,n introduced in 



Definition E31 

Furthermore, we show that the fact that a function u belongs to the space D(A k ) can be 
formulated in terms of the decay rate of Vu. These results will be useful for the resolution of 
evolution problems involving the spatial operator A as well as for the analysis of the properties 
of a solution. For this part, we follow Section 4 of [2\. 

Definition 7.1. i) Let / C K be a closed interval and let a : / — > M be a measurable, 
non-negative function with <r(A) > for all A G /. Define L 2 (I, a) by 



(F,G)v M := l a(X)F(X)G(X)dX, 

r) 



\F\V>{I,*) ■— (F, F ) l L?{I,a) 



L 2 (I,a) := {F : I — ► C measurable, \F\ L 2r I>(7 -\ < +oo}. 
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ii) Consider now the weights a k '■= Qk f° r k = 1, . . . ,n, where q k is given in Theorem 15.31 
We endow l? a := IJfc=i L 2 {[ a k-, +o°), with the inner product 

n 

(F, G) a := ^(i ? fc,G , fc ) i 2 ([afei+oo)i(J;c) 

fc=l 

and denote 

\F\. := (F,F) l J 2 . 

Note that I? a is a Hilbert space, since it is the product of the Hilbert spaces L 2 ([a k , +oo), a k ). 

Now, we define the transformation V, together with its (right) inverse Z, which is later on 
shown to diagonalize A. 

Definition 7.2. i) For all / G L 1 (A r , C) and k = 1, . . . , n the function (Vf) k : [a k , +oo) — ► 
C is defined by 



(Vf) k (X) := / f(x)(F-^)(x) dx. 



N 



ii) Consider \ £ G°°(R) such that x = on (— oo, a n + 1) and x = 1 on (a n + 2, +oo). For 
Gfc € C°°([afc, +oo), C), such that ^Gfc can be extended by zero to an element of <S(R), 
k G {1, . . . , n}, we define Z(G X , . . . , G n ) : N -> C by 

. . . , G n )(x) == X) / ^(A)G fc (A)(F-' fc )(x) dA. 

k=1 J (a k ,+oo) 

Note that in contrast to the considerations in [2], the functions (Vf) k here are complex- valued, 
due to the choice of the generalized eigenfunctions. 

Lemma 7.3. Consider f G L 2 (N) with compact support, x us in Definition |7.£| and for k G 
{1, ...,n} let again G k G G°°([a fc , +oo)) fre suc/i i/iai xGk G <S(M). TTien G = (Gi, . . . ,G n ) G 
L 2 , G L\, Z{G) G L 2 (iV) = H and 

(G,Vf) CT = (Z(G)J) H . 

Proof. We have 

n n „ 

(G,Vf) a = J2{G k , (Vf) k ) L2 ,, + , k) = / a k (X)G k (X)(VfTkW dX 



fc=l fc=l 



V / a fc (A)G fc (A)( / /(x)(F"' fc )(x) dx) dA 




fc=1 J(a k ,+co) 

a k (X)G k (X)(F~' k )(x) dXjf{x)dx = (Z(G),f) H . (16) 

It only remains to make sure that the assumptions for Fubini's Theorem are satisfied to justify 
(fT6|) . In fact, it is sufficient to estimate (A,x) i— > a k (X)(F^ ,k )(x) on [ajt,ajfc+i] X N for a fixed 
k G {1, . . . , n}, since / is compactly supported and G k is rapidly decreasing. 

In order to do so, recall that c k £ k (X) = \/ c k (X — a k ) belongs to R if and only if A > a k and 
to iR otherwise and that for any A > a%, we have |u>(A)| 2 > J2?=i c i\^ ~ a i\ due to Lemma [3771 



N k=1 J(a k ,+oo) 
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Thus, putting in the expression for Sfc(A), cf. Definition 13.41 we get for x £ N k 



\a k (X)(F-' k )(x)\ 



cos 



< 



\w(X)\ 2 
|cfc6(A)| 



(Cfc(A)a?) - is k (\) sin(£ k (X)x)) 



< 



Er=i c /I A - a ;l 



:i+k(A)D 



Er=iQi6(A)i < zr=i 



Er=i c i|A-a;| |c fc ^fc(A)| 
Furthermore, for x £ Nj, j > k, we have 



MA)(F-< fe )(*)| 

Finally, for x £ Nj, j < k, 

\a k (\)(F-> k )(x)\ 



w(A)| s 



CfcCfc(A) (A)a: 



< 



Er=i c dA-a«| 



cxp 



Cfc\/|A - a*. 



Cfc4(A) 



|to(A)|s 



< 



For all three cases the bound is a continuous function of A on [afc] a fc+i] and so, belongs to 
Lf oc ([a k ,a k+ i]), which is enough for Fubini's Theorem, due to the properties of / and G k . □ 

Lemma 7.4. Consider f £ Yl k =i Cc°(-^fc)- Then, for every k £ {1, . . . , n}, 

(Vf) k £ C([a k , +oo)) n C°°([a fc , a k+1 )) n . . . D C°°(K, +oo)) 

anc? X°~k(Vf)k £ <S(R) with x as in Definition \ 7. 2\ 



Proof. As in Section 4 of [2], (Vf) k is a linear combination of Fourier and Laplace transforms 
of functions in C%°(Nj), j £ {1,... , n}. □ 

The next step is to show that V is an isometry and Z is its right inverse. The proof of the 
following theorem is precisely as in [2J. 

Theorem 7.5. Endow Wl =l C^{N k ) with the norm of H = \\l =l L 2 {N k ). Then 

i) y : Ilfc=i Cc°(-^fc) — * L 2 is isometric and can therefore be extended to an isometry 
V:H^L 2 a . 

In particular, for all f £ H we have \Vf\ 2 = (/, /)h- 

ii) Z = V- 1 on V(H). 

iii) Z can be extended to a continuous operator on L 2 . 

Note that it is not clear, whether Z is injective, in contrast to the situation for n = 2, see [2j. 
There the proof of injectivity relies on a choice of generalized eigenfunctions that are real- valued. 
Here, a choice of such eigenfunctions would probably destroy the feature that the matrix q, found 
in Theorem 15.31 is diagonal. And it is this property that will allow us to carry over from [2J the 
proof of Theorem 17.71 However, the injectivity of Z is not important for the applications we 
have in mind. 

In the sequel we shall again write V for V for simplicity. 

Our final aim is to show that V diagonalizes the operator, i.e. V \A J 'u)(A) = A J V(n)(A) for all 
u £ D{A :> ). In order to formulate this precisely for a measurable function tp : R — > R we denote 
the corresponding multiplication operator on L 2 by M^, i.e. for a function F £ L 2 we have 

(M^F) fe (A) = V(A)F fc (A), AG[a fc ,+oo), k £ {1, . . . , n}. 

Then we have the following lemma, whose proof is again analogous to [2J Lemma 4.12]. 
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Lemma 7.6. i) Let j £ N and pj : R — > R 6e defined by Pj(x) = x 3 . Then for any 
f £ -D(yP') we have 

V(A 3 f) = M Pj Vf. 

ii) Lei Vl/ : [ai,oo) -> 1 ie a bounded, measurable function. Then *&(A) defined by the 
spectral Theorem is a bounded operator on H and for all f £ H we have 

V(V{A)f) = M 9 (Vf). 

Finally, this leads to a characterization of the spaces D(A 3 ). 

Theorem 7.7. For j £ N the following statements are equivalent: 

i) u € D(AP), 

ii) \»X>(Vu)(\) eL 2 a , 

iii) A i-» X j (Vu) k (X) £ Lj([ofc, +oo)), fc = l,...,n. 

Proof. Due to Theorem 17. 5 [ it holds (A 3 u, A 3 u)h = \V(A J u)\^. Now Lemma 17.61 implies 

(A>u,A j u) H = \M Pj Vu\l. 

□ 

The above results provide explicit solution formulae for evolution equations involving the 
operator A. For example for the wave equation u + Au = with u(0) = uq and u(0) = 
a formal solution is given by u(t) = Z cos(v / At)V r no. Our aim in the future will be to study 
properties of these solutions, as indicated in the introduction. 

8. Appendix 

The book of J. Weidmann [27] describes a general approach to the spectral theory of systems of 
Sturm-Liouville equations, which is in principle applicable to our setting. In this point of view, 
our problem is seen as a system of n equations on (0, +oo) coupled by boundary conditions in 
0. Thus the kernel of the resolvent is a matrix-valued function JC(-, •, A) : (0, +oo) x (0, +oo) — > 
C nxra . The relation to Definition 13.61 is given by 

ICij(x,x',X) = K(x,x',X) where (x,x') £ Ni x Nj = (0, +oo) x (0, +oo). 

The fundamental hypothesis of Weidmann is that the kernel of the resolvent can be written in 
the following form: 

p 

J2q=l [S2 a lq W l( X > A )) WqW, A) T , for x' < X, 



JC(x, x , A) 



l=i 



V 

^?=1 (X)^ W l( X ' X i) w q( x 'i X ) T i for x> > x > 
1=1 



(17) 



=: m,q{x,\) 

where ai q ,/3i q £ C and the w q : [0, oo) x C — > C n are such that {w q (-,X) : q = l,...,p} 
is a fundamental system of ker(Af — XI), the space of generalized eigenfunctions of A. Here 
Af : D(Af) — ► C (R) is the formal operator, in our case Af = A = (— • d% + ak)k=i,...,n but 
D(Af) = YYk = i C 2 (Nk) , i.e. the operator A without transmission conditions nor integrability 
conditions at oo. Clearly in our case dim(ker(Aj — A/)) = 2n and thus p = 2n. 

In contrast to the typical applications treated in [27J as for instance the Dirac system, we 
consider in this paper a transmission problem, i.e. intuitively the components of the w q are 
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functions on different domains N k (while mathematically all N k are equivalent to (0,+oo)): 
the branches of a star. For all such applications it is highly important to use only generalized 
eigenfunctions satisfying the transmission conditions (To) and (T\), for example Theorem 17.71 
would be impossible otherwise. 

Supposing the ansatz (|17p of Weidmann, this is not possible, what we shall show now. 

In Definition 13. 61 we have given an explicit expression for the (unique) kernel K of the resolvent 
of A, using only elements of kei(AT-XI), where At : D(At) — ► H satisfies At = A = (— cu-d 2 + 
a fc )fc=l,...,n and D(A T ) = JILi C 2 {N k )C\{{u k )l =1 satisfies (T ), (Ti)}. Note that kei{A T - XI) is 
the n-dimensional space of generalized eigenfunctions of A satisfying the transmission conditions 
(To) and (Ti), but without integrability conditions at infinity. 

Suppose that we have a representation of K as given in (fTT|) satisfying 

w q (; A) G ker(^ r - XI), q=l,...,p, AG p(A) 

and thus we can take p = n. Let us fix x G N\ and A G p(A). Then the rriq(x, X) =: rriq G C are 
constants and the expression 

v 

g(x', X) := ^2 mqW q (x', X) 

9=1 

defines a function 

g(-,X) : {x' G [0,oo) : x > x} -»■ C n . 

Clearly, the functions 

[«/,(•, A)k:J\^-C 

can be uniquely extended to entire functions in x' because they are linear combinations of 
e ±*€jW x \ Thus the [g(-, X)]j are entire. Our hypothesis w q (-, X) G ker(AT — XI), q = 1, . . . ,p, 
implies thus 

g(-,X) eker (A T -XI), q = l,...,p. (18) 
But comparing (|17|) with Definition 13.61 yields 

g{x', X) = F±' 2 (x) • F±'\x') for x' G A^ 2 , . . . , N n 

and 

g{x, X) = F^' 1 (x) ■ F^ 2 (x) for x' G N\, 

which makes sense after analytic continuation. Using the assumption that g satisfies (To), we 
get from these two equalities, putting x' = G A^2 into the first and x' = G N% into the second 

F±> 2 (x) = F^ 2 (x)F^(0) = 5 (0, A) = F^(x)F^ 2 (0) = F^(x). 

But inspecting the definitions of the generalized eigenfunctions, bearing in mind that x G Ni 
was arbitrary, this implies 

cos(£i(A)x) ± i sin(£;i(A)x) = cos(£i(A)x) ± isi(X) sin(£i(A)x) 

and thus si(A) = 1. This finally implies Y^=i c j£jW = 0' which is impossible for all real A < Oi, 
since then £j(A) is purely imaginary. Furthermore, Y^j=i c j£j ^ s analytic, so this equality can, if 
ever, only be fulfilled on a discrete set of A G C. 
We have thus proven: 

Theorem 8.1. Representation (|17|) of the kernel JC of the resolvent of A is not possible, using 
exclusively generalized eigenfunctions 

w g (-, X) G kei(A T -XI), q = l,...,p. 
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The reason for this is, shortly speaking, the rigidity of (|17p caused by the use of the same 
linear combinations of generalized eigenfunctions above the diagonals of all Nj x N/. (idem 
below). This is not compatible with the fact that the kernel is non-smooth only on the main 
diagonals, i.e. the diagonals of Nj x Nj, cf. Definition 13.61 and Figure [TJ 

This means that the approach of J. Weidmann in [27], when applied to problems on the star- 
shaped domain, does not sufficiently take into account its non-manifold character. The resulting 
expansion formulae would use generalized eigenfunctions which are in a sense incompatible with 
the geometry of the domain. This would have undesirable consequences: for example, the 
important feature of Theorem 17.71 that the belonging of u to D{A 3 ) is expressed by the decay 
of the components of Vu would be impossible, due to artificial singularities in the expansion 
formula. 
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